Aging dynamics of heterogeneous spin models 

Andrea Montanari 

Laboratoire de Physique Theorique de VEcole Normale Superieure* 
24, rue Lhomond, 15231 Paris CEDEX 05, FRANCE 
Internet: Andrea. MontanariOlpt .ens .fr 

Federico Ricci-Tersenghi 

Dipartimento di Fisica and SMC and UdRl of INFM 
Universita di Roma "La Sapienza" 
Piazzale Aldo Moro 2, 1-00185 Roma, ITALY 
Internet: Federico . RicciOromal . inf n . it 



May 1, 2003 



Abstract 

We investigate numerically the dynamics of three different spin models in the aging regime. 
Each of these models is meant to be representative of a distinct class of aging behavior: 
coarsening systems, discontinuous spin glasses, and continuous spin glasses. In order to study 
dynamic heterogeneities induced by quenched disorder, we consider single-spin observables for 
a given disorder realization. In some simple cases we are able to provide analytical predictions 
for single-spin response and correlation functions. 

The results strongly depend upon the model considered. It turns out that, by comparing 
the slow evolution of a few different degrees of freedom, one can distinguish between different 
dynamic classes. As a conclusion we present the general properties which can be induced from 
our results, and discuss their relation with thermometric arguments. 
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1 Introduction 



Physical systems with an extremely slow relaxation dynamics (aging) are, at the same time, 
ubiquitous and fascinating [Q. Much of the insight we have on such systems comes from the 
study of mean-field models |2j. 

One of the weak points of the results obtained so far is that they focus on global quantities, 
e.g. the correlation and response functions averaged over all the spins. On the other hand, we 
expect one of the peculiar features of glassy dynamics to be its heterogeneity PlilloTloTITllSllMTUI 
lllj . In order to understand this character, we study the out-of-equilibrium dynamics of three 
models belonging to three different families of slowly evolving systems: coarsening systems, 
discontinuous and continuous glasses. 

The correlation and response functions of a particular spin depend (among the other things) 
upon its local environment, i.e. the strength of its interactions with other spins. However, 
the way the single-spin dynamics is influenced by its environment is highly dependent upon 
the nature of the system as a whole. For instance, as we will show, while in coarsening 
systems strongly interacting spins relax faster, for discontinuous glasses the opposite happens. 
Continuous glasses lie somehow midway. In principle, this allows to distinguish different types 
of slow dynamics just by looking at the relation between a couple of spins. 

The local environment is not the only source of heterogeneities |12U13lll4j . Even in systems 
without quenched disorder, like structural glasses, the thermal noise is able to break the initial 
spatial uniformity and to bring the system in a strongly heterogeneous configuration [HHUO)]. 
Nevertheless we think that considering systems with quenched disorder can be an instructive 
first step even in that direction. As it has been argued several times |15|llri| . structural glasses 
behaves similarly to some disordered systems because of a sort of self-induced disorder. Each 
molecule relaxes in the amorphous environment produced by the (partially frozen) arrangement 
of the other ones. 

In order to extract quenched-disorder-induced heterogeneities we will average on a very 
large number of independent thermal histories. This will delete the effects of the thermal 
noise. 

On the contrary we are forced not to perform a naive average over the disorder realizations, 
because this would wash any difference between spins. Instead of doing more careful disorder 
averages, e.g. conditioning on the local environment of the spin under study, we prefer to 
work with a unique fixed disorder realization. In the limit of large system size we expect local 
quantities still to fluctuate from site to site, and to converge in distribution sense, making the 
analysis of a single typical sample representative of the whole ensemble. 

After these preliminaries, we can summarize the approach used in this work. Given a 
disordered model, we take a few typical samples 1 (as big as possible according to our numerical 
capabilities) from the ensemble and we repeat a huge number of times the typical numerical 
experiment used for studying out-of-equilibrium dynamics: start from a random configuration, 
quench the system to a low temperature, where it evolves slowly, wait a time t w , switch 
on a small perturbing field, and take measurements. The observables we measure are local 
quantities, like single-spin correlation and response functions, averaged over the thermal noise. 

We shall consider three different disordered models: 

• a two-dimensional ferromagnetic Ising model (couplings are all ferromagnetic but of dif- 
ferent strengths), which has a ferromagnetic phase below the critical temperature; 



1 The complete specification of the samples used in our simulations is available, upon request to one of the authors. 
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• the 3-spin Ising model on random hypergraph, which has a glassy phase with one step of 
replica symmetry breaking (1RSB); 

• the spin glass Ising model on random graph, also known as Viana-Bray (VB) model |17j . 
which is believed to have a glassy phase with continuous replica symmetry breaking 
(FRSB). 

The last two models are examples of diluted mean-field spin glasses. They lack of any 
finite-dimensional geometric structure: this makes them soluble using mean-field techniques. 
On the other hand, the local fluctuations of quenched disorder are not averaged out as in com- 
pletely connected models. For instance, the local connectivity is a Poissonian random variable. 
Because of these two features, they are an interesting playing ground for understanding het- 
erogeneous dynamics. 

Diluted mean-field models have been intensively studied in the last years, one of the qualify- 
ing motivations being their correspondence with random combinatorial problems ^B]- Statical 
heterogeneities have been well understood, at least at 1RSB level. Throughout the paper 
we shall neglect FRSB effects ^2], and assume that 1RSB is a good approximation. In 
Refs. |2(JU21| . the authors defined a linear-time algorithm that computes single-spin static 
quantities for a given sample in 1RSB approximation. The algorithm was dubbed surveys 
propagation (SP) and, strictly speaking, was defined for computing zero-temperature quan- 
tities. It is straightforward, although computationally more demanding, to generalize it for 
finite temperatures T (the generalization follows the ideas of Ref. H2]) : we shall call this 
generalization SPy. 

The resulting heterogeneities can be characterized by a local Edwards- Anderson parameter. 
This can be defined j^Sl by considering m weakly-coupled "clones" {c^ 1 ), . . . , <j( m )} of the 
system. The local overlap between two of them q^l_(m) = (<r t af ), with a ^ b, is given by 

^(-)-^ E e-^-MS, (1-1) 

a£ {states} 

where the sum on a runs over the pure states, (-) Q denotes the thermal average over one of 
such states, and Z m = ^ a e~^ mFa . Equation follows from the observation that the m 
clones stay at any time in the same state a, and that each state is selected with probability 
e~P mFa / Z m . The parameter m enables us to select metastable states. In fact we expect 
the dynamics of discontinuous glasses to be tightly related with the properties of high-energy 
metastable states (SUES]- 

While in a paramagnetic phase (/ea(^) = apart for a non-extensive subset of the spins, 
in the spin glass phase (/EA( m ) > in a finite fraction of the system. In general q$ A {m) 
depends upon the site i: the phase is heterogeneous. We will return in the next Section on 
the dynamical significance of this and other statical results. 

The paper is organized as follows. In Sec. Elwe present some of the theoretical expectations 
which we are going to test. We also give a few technical details concerning the numerics. 
Section deals with coarsening systems. We postulate the general behavior of response and 
correlation functions, and test our predictions on a simple model. In Sees. 0] and |S] we 
present our numerical results for, respectively, the 3-spin and 2-spin interaction spin glasses 
on random (hyper) graphs. The particularly easy case of weakly interacting spins is treated in 
Sec. El We show that the aging behavior of these spins can be computed from the behavior 
of their neighbors. Finally, in Sec. [7J we discuss the general picture which emerges from our 
observations. In Sec. [HI we interpret some of these properties using thermometric arguments. 
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Appendix ^present some calculations for coarsening dynamics. A brief account of our results 
has appeared in Ref. 



2 Generalities 

In the following we shall discuss three different spin models. Before embarking in such a tour 
it is worth presenting the general frame and fixing some notations. 

Our principal tools will be the single-spin correlation and response functions: 



Cij(t,t w ) = (<Ji{t)(Tj{t w )) , Rij(t,t w ) 



dhj(t v 



(2-1) 

h=0 



where the average is taken with respect to some stochastic dynamics, and hj is a magnetic 
field coupled to the spin j. It is also useful to define the integrated response Xij{t,tw) = 

We shall not repeat the subscripts when considering the diagonal elements of the above func- 
tions (i.e. we shall write Cj for Cu, etc.). The global (self-averaging) correlation and response 
functions are obtained from the single-site quantities as follows: C(t, t w ) = (1/N) ]T\ d(t, t w ); 
x(t,t w ) = (l/N)Y^iXi(titw)- The times t and t w are measured with respect to the initial 
quench (at £ que nch = 0) from infinite temperature. 

We will be interested in comparing the outcome of static calculations and out-of-equilibrium 
numerical simulations. For instance, we expect the order parameter to have the following 
dynamical meaning 2 

?BA(rath) = lim lim Ci(t w + At, t w ) , (2.2) 

At— >oo t w ^QO 

where m t h is the parameter which select the highest energy metastable states. 

In the aging regime At,t w ^> 1, Ci(t w + At,t w ) < (^(mth)- We expect the functions (|2.1|) 
to satisfy the out-of-equilibrium fluctuation-dissipation relation (OFDR) [23] 



TR^t,^) = ^[CiCMuOl^CiCMuO • (2-3) 

If -Xj[Cj] = 1 the fluctuation-dissipation theorem (FDT) is recovered. The arguments of 
Refs. [2*7|l28| . and the analogy with exactly soluble models [2*Hl29p;-i()j suggest that the function 
Aj[C] is related to the static overlap probability distribution: 

w = -^. (-) 

For discontinuous glasses the dynamics never approaches thermodynamically dominant states. 
In this case the function Pi(q) entering in Eq. (|2.4|) is the overlap distribution among highest 
metastable states. We refer to the next Sections for concrete examples of the general relation 

Let us now give some details concerning our numerical simulations. We shall consider 
systems defined on iV Ising spins crj = ±1, i S {1, ...,N}, with Hamiltonian H(a). The 
dynamics is defined by single spin flip moves with Metropolis acceptance rule. The update 



2 Notice that we shall always work with finite samples. The limits At, t w — > oo must therefore be understood as 
1 <C At, t w <C t crg (N). After a time of order t e i S (N) ergodicity is re-established and the system equilibrates. 



4 



will be sequential for the spin glass models of Sees. |1] and |S] and random sequential for the 
ferromagnetic model of Sec. |3J 

For each one of the mentioned models, we shall repeat the typical aging "experiment". The 
system is initialized in a random (infinite temperature) configuration. At time t qU ench = 0, the 
system is cooled at temperature T within its low temperature phase. We run the dynamics 
for a "physical" time t w (corresponding to t w attempts to flip each spin). Then we "turn on" 
a small random magnetic field hi = ±ho and go on running the Metropolis algorithm for a 
maximum physical time AiMAX- Notice that the random external field is changed at each 
trajectory. 

The correlation and response of the single degrees of freedom are extracted by measuring 
the following observables: 

. t u ,+2At 

Ci(t w + 2At,t w \ha) = — Yl ( 2 - 5 ) 

t'=t w +At+l 
. t w +2At 

Xi (t w + 2At,t w \h ) = Yl (vi(t')sim(hi)), (2.6) 

t'=t w +At+l 

where (■) denotes the average over the Metropolis trajectories and the random external field. 
The sum over t! has been introduced for reducing the statistical errors. While it is a drastic 
modification of the definition (|2.1j) in the quasi-equilibrium regime At <C t w , it produces just 
a small correction in the aging regime At,t w 3> 1. This correction should cancel out in two 
interesting cases: (i) in the time sector t/t w = const., if one restrict himself to the response- 
versus-correlation relation; {ii) in "slower" time sectors (e.g. (logi)^ — (\ogt w )^ = const., 
with fi < 1). The functions (|2.5j) and (|2.6j) have finite ho — > limits Ci(t w + At,t w ) and 
Xi{t w + At,t w ). 

Finally, let us mention that we shall look at the Xi(t^w) versus Ci(t,t w ) data from two 
different perspectives. In the first one we focus on a fixed site i and vary the times t and t w : 
this allows to verify the relations ()2.3|) and (|2.4|) . We shall refer to this type of presentations 
as FD plots. In the second approach we plot, for a given couple of times, all the points 
(Ci(t, t w ),Xi{t-> tw)) for i = 1, . . . , N. Then we let t grow as t w is kept fixed. We dubbed such 
a procedure a movie plot. It emphasizes the relations between different degrees of freedom in 
the system. 

3 Coarsening systems 

Coarsening is the simplest type of aging dynamics |32U31| . Despite its simplicity it has many 
representatives: ferromagnets (both homogeneous and not), binary liquids, and, according to 
the droplet model [3*3*ll3Hl3l)ll3l)] . spin glasses. 

Consider a homogeneous spin model with a low temperature ferromagnetic phase (e.g. an 
Ising model in dimension d > 2). When cooled below its critical temperature, the system 
quickly separates into domains of different magnetization. Within each domain the system is 
"near" one of its equilibrium pure phases. Nevertheless it keeps evolving at all times due the 
the growth of the domain size £(t). This process is mainly driven by the energetics of domain 
boundaries. 

In the t — > oo limit, the coarsening length obeys a power law £(i) ~ t l l z (for non-conserved 
scalar order parameter z = 2). Two-times observables decompose in a quasi-equilibrium part 
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describing the fluctuations within a domain (C eq and Xeq in the equations below), plus an 
aging contribution which involves the motion of the domain walls (C ag , Cd w and Xdw) : 

C(t,t w ) ~ C eq (i — t w ) + q EA Ca,g ( ^r: \ ) ~ 9EAi M , a Cd w ( j\ \) , (3-1) 

x(t,t w ) « xcq (*-*«,) + + W^Xdw (^Py) » ( 3 - 2 ) 

where (jea is the equilibrium Edwards- Anderson parameter. For a ferromagnet qEA = M((3) 2 , 
M(P) being the spontaneous magnetization. Moreover C eq (r) decreases from (1 — <?ea) to 
as r goes from to oo, and C ag (A) goes from 1 to as its argument increases from 1 to oo. 
Finally the equilibrium part of the susceptibility Xeq(i~) goes from to (1 — qsp^/T. In the 
case of a scalar order parameter both the response and correlation functions receive subleading 
contributions (Cdw and Xdw) from spins "close" to the domain walls. Notice that these spins 
will decorrelate faster and respond easier than the others (in other words Cd w and Xdw are 
typically positive). These contributions are expected to be proportional to the density of 
domain walls Pd w (t w ) (x £(t w ) . This would imply o = a' = 1/z. 

It is easy to generalize this well-established scenario to include single-spin quantities in 
heterogeneous systems. We expect that the quasi-equilibrium parts of the correlation and 
response functions will depend upon the detailed environment of each spin. On the other 
hand, the large scale motion of the domain boundaries will not depend upon the precise point 
of the system we are looking at. Therefore the aging contribution will depend upon the site 
x only through the local Edwards- Anderson parameter q x A = M x (f5) 2 . The reason is that q x A 
quantifies the distance between pure phases as seen through the spin a x . 

We are led to propose the following form for single-site functions: 

C x (t,t w ) « C x q (t — t w ) + q x A C ag ^yl — Q^: A t w a Cdw ^jl ' 
Xx (t, t w ) « X ?(t ~t w )+ + q* A t- a ' Xdw (4¥K) ■ (3.4) 



This ansatz can be summarized, as far as 0(t~ a ,t~ a ') terms are neglected, in the schematic 
response-versus-correlation plot reported in Fig. ^ Each spin follows its own fluctuation- 
dissipation curve. This is composed by a quasi-equilibrium sector T X x = 1 — C x , plus an 
horizontal aging sector T\x = 1 — q x A - Moreover, for each couple of times t w and i, all the 
points are aligned on the line passing through (C = 0,Tx = 1). 

Notice that spins with larger q x A relax faster. Roughly speaking this happens because they 
have to move a larger distance in order to jump from one pure phase to the other. 

Due to the independence of Cd w and Xdw upon the site x, Eqs. (|3.3|) and (|3.4j) imply that 
alignment in the x~C plane is verified even in the pre-asymptotic regime. This property is 
therefore more robust than the OFDR which is violated by 0(t~ a ,t~ a ) terms, cf. Figs. 0]and 

El 

Finally, both the large-n calculation of Sec. 13.1.11 and our numerical data, cf. Sec. 13.1.21 
suggest that domain- wall contributions in Eqs. (|3.3|) . ()3.4|) have the same order of magnitude 
Cdw (A) ~T Xd w(A). 
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Figure 1: Qualitative picture of the response- versus-correlation plot for coarsening systems. Bold 
dashed lines are the single-spin OFDR's. Black dots represent the correlation and response functions 
at a given pairs of times (t,t w ). Arrows correspond to the "velocities" of the dots when they move 
along the fluctuation-dissipation curves. 



3.1 A staggered spin model 

Here we want to test our predictions in a simple context. We shall consider a d-dimensional 
lattices ferromagnet, defined by the Hamiltonian 

where the sum runs over all the couples (xy) of nearest-neighbors on the lattice and J xy > 0. 
Moreover we assume periodicity in the couplings. Namely, there exist positive integers li,...,ld 
such that, for any x, y G and € {1, . . . , d} 

Jxy = Jx+p,ln,y+tdf, (3-6) 

where p, is the unit vector in the fi-th direction. Clearly there are V = l\ ■ I2 • •• Id different 
"types" of spins in this model. Two spins of the same type have the same correlation and 
response functions. We can identify these V types with the spins of the "elementary cell" 
A = {x G Z d \0 <x^< l^}. 

Spatial periodicity is helpful for two reasons: (i) it allows an analytical treatment in the 
large-n limit; (ii) averaging the single-spin quantities over the set of spins of a given type 
greatly improves the statistics of numerical simulations. 

3.1.1 Large n 

The model (|3.5|) is easily generalized to n- vector spins cj> x = (<%., ...,<$£). We just replace the 
ordinary product between spins in Eq. (|3.5|) with the scalar product. Moreover we fix the spin 
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Figure 2: Domain wall contributions to the correlation (dashed line) and integrated response func- 
tions (continuous line) in the n — > oo limit: the universal scaling functions, cf. Eqs. (j3.11|) and 
(j3.12J) . as a function of A = ^Jt/t w . 



length: <p x ■ <j> x = n. The dynamics is specified by the Langevin equation: 

dt€(t) = -Ut)€(t) + E J ^<W + vZ® , (3-7) 

y 

where we introduced the Lagrange multipliers Cx(t) in order to enforce the spherical constraint. 
The thermal noise is Gaussian with covariance 

(Vx(t)v b y (s)) = 2T5 xy 5 ab 5(t - s) . (3.8) 

The definition of correlation and response functions must be slightly modified for an n- 
components order parameter: 

Like its homogeneous relative |31| . this model can be solved in the limit n — > oo. The 
calculations are outlined in App. Let's summarize here the main results. For d > 2 the 
model undergoes a phase transition at a finite temperature T c . Below the critical temperature 
the 0{n) symmetry is broken: ((/>".) = M x ((3)5 al . Of course the spontaneous magnetizations 
preserves the spatial periodicity of the model: M x +i = M x {(3). 

At low temperature, the forms (j33|) - (j3~l)) hold, with q^ A = M%(J3), a = a' = d/2 - 1, z = 2 
and 



Cag(A) = [—^- ) ■ (3-10) 
The subleading contribution read 

2T 

° dM) = p^i^'^ 4 (3 - 11} 

Xdw(A) = (47r)rf/2(E J gAM 2 )Al/2 ^x(A;rf), (3.12) 
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Figure 3: Definition of the ferromagnetic couplings for the two-dimensional model studied in 
Sec.EU 



where A is a constant which depends uniquely on the couplings J xy , cf. App. IA.1I ^ r c(-), 
J~ x {-) are two universal function which do not depend either on the temperature or on the 
particular model. The explicit expressions for these functions are not very illuminating. We 
report them in App. IA.21 see Eqs. (|A.18|) and ()A.19|) . Here we plot the two functions in the 
d = 3 case, see Fig. |2I Notice that both .Fc(A;(f) and T x (\;d) vanish in the A — > oo limit. 
This could be expected because we know that C x (t,t w ) — > and Xx(t,t w ) — > (1 — Mj)/T as 
t — > oo for any fixed t w . 

3.1.2 Numerical simulations 

We simulated the model (|3.5[) in d = 2 dimensions with l\ = li = 2 and the choice of couplings 
among spins in the elementary cell illustrated in Fig. 13 We used square lattices with linear 
size L. There are V = 2 2 different type of spins in this case. We shall improve our numerical 
estimates by averaging the single site functions Ci(t,t w ; ho) and Xi{t->tw\ ho), cf. Eqs. 1)2. 5 Jl 
and (|2.5|) . over the L 2 /4 spins of the same type. 

Most our numerical results were obtained at temperature T = 1. A rough numerical 
estimate yields T c = 1.10(5) for the critical temperature. The equilibrium magnetizations for 
T = 1 of the four type of sites are: M = 0.8803(5), M l = 0.8395(5), M 2 = 0.7573(5), and 
M3 = 0.8624(5). Notice that, in order to separate the magnetization values on different sites, 
we are forced to choose a quite high temperature for our simulations. 

We expect the growth of the domain size in the model (|3.5j) to follow asymptotically the 
law £(i) k(f3) ■ t x / z , with z = 2, as in the homogeneous case. The pinning effect due 
to inhomogeneous couplings will renormalize the coefficient k((3). We checked this law by 
studying the evolution of the total magnetization starting from a random initial condition for 
different lattice sizes. It turns out that the law is reasonably well verified with a coefficient 
k((3 = 1) of order one. 

The aging "experiment" was repeated for several values of the waiting time t w = 10, 10 2 , 
10 3 , 10 4 , 10 5 . The correlation and response functions were measured up to a maximum time 
interval (respectively) A^max = 2 10 , 2 13 , 2 15 , 2 17 , 2 19 . The linear size of the lattice was 
L = 2000 in all the cases except for t w = 10 5 . In this case we used L = 1000. All the results 
were therefore obtained in the £(i) <C L regime, with the exception, possibly, of the latest 
times in the t w = 10 5 run. Some systematic discrepancies can be indeed noticed for these 
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Figure 4: Correlation function and FD plot for type-0 sites (cf. Fig. |HJ). Different symbols corre- 
spond to t w = 10 (□), 10 2 (*), 10 3 (x), 10 4 (+), 10 5 (filled □). The dot-dashed line on the left is the 
equilibrium Edwards- Anderson parameter Mg. On the right we report the FDT line T\ = 1 — C 
(dot-dashed) and the OFDR (dashed) which corresponds to Eqs. ()3.3|) and (|3.4|) . 



data. In the Table below we report the number iV s t a t of different runs for each choice of the 
parameters. 
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Let us start by illustrating how the asymptotic behavior summarized in Fig. Q is ap- 
proached. In Fig. 0]we show the correlation functions and the FD plot for type-0 sites. Notice 
that the approach to the asymptotic behavior is quite slow and, in particular, the domain-wall 
contribution to the response function is pretty large. This can be an effect of the proximity 
of the critical temperature: the "thickness" of the domain walls grows with the equilibrium 
correlation length. 

In Fig. Owe verify the alignment of different sites correlation and response functions for a 
given pair of times (t,t w ). Notice that the alignment works quite well even for "pre-asymptotic" 
times, i.e. when the anomalous response is still sizeable and the OFDR is not well verified, cf. 
Fig. H 

In order to check the form (|3.4|) for the site dependence of the domain wall contribution, 
we plot in Fig. El the rescaled response and correlation functions: 

d eS = ^C x , T X T = l--^(l-T Xx ), (3.13) 

Six Six 

where q is an arbitrary reference overlap. The rescaled correlation and response functions of 
different types of spin coincide perfectly for any couple of times (t,t w ). 

Finally, we notice that we can consistently define a time-dependent fitting temperature as 
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Figure 5: Movie plots at t w = 10 3 (left) and 10 4 (right). The various symbols correspond to different 
types of spin: type (+), type 1 (x), type 2 (□), and type 3 (*). The straight lines confirm the 
alignment predicted in the general picture, cf. Fig. [T] 
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Figure 6: Correlation function (left) and FD plot (right) with the rescaled correlation and response 
functions, see Eq. f!3.13|) . for all the four spin types and several different waiting times: t w = 10 2 , 
10 3 , 10 4 , and, on the left 10 5 . Here h = 0.05. 
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Figure 7: The fitting temperature (j3.14|) as a function of A = ^/t/t w for t w = 10 (filled □), 10 2 
(□), 10 3 (*), 10 4 (+), and 10 5 (x). The dot-dashed line is the n = oo scaling function (ETTU1) . with 
d = 2. 



the slope of the lines in Fig. [""J i.e.: 

T*(t,t w )= T °* (t >^ (3.14) 
1 - Txx(t,t w ) 

As a consequence of Eqs. ()3.3[1 and (|3.4j) this temperature should depend upon t and t w only- 
through the parameter A = £(t)/£(t w ). In Fig. [JJwe verify this scaling. 

4 Discontinuous glasses 

In this Section we consider a ferromagnetic Ising model with 3-spin interactions, defined on a 
random hypergraph [3*7113*5] . More precisely, the Hamiltonian reads 

(ijk)en 

The hypergraph Ti defines which triplets of spins do interact. We construct it by randomly 
choosing M among the N(N — 1)(N — 2)/3! possible triplets of spins. 

Although ferromagnetic, this model is thought to have a glassy behavior, due to self- 
induced frustration [39] . Depending upon the value of 7 = M/N, it undergoes no phase 
transition (if 7 < 7^), a purely dynamic phase transition (if 7^ < 7 < 7 C ), or a dynamic and 
a static phase transitions (if 7 > 7 C ) as the temperature is lowered. The 1RSB analysis of 
Refs. [H3IIII! yields 7^ w 0.818 and 7 C sa 0.918. These results have been later confirmed by 
rigorous derivations [4*l*]l4*2]. 

We studied two samples extracted from the ensemble defined above: the first one involves 
N = 100 sites and M = 100 interactions (hereafter we shall refer to it as Ha)', in the second 
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Figure 8: The complexity E(T) (left) and the 1RSB parameter m^iT) (right) for threshold states 
as functions of the temperature T. These curves refer to sample Ha considered in Sec. HI 



one (Hb) we have N = M = 1000. In both cases 7 = 1 > j c . The hypergraph Ha consists 
of a large connected component including 96 sites, plus 4 isolated sites (namely the sites 
i = 15,22,62,69). The largest connected component of Hb includes 938 sites (there are 62 
isolated sites). We will illustrate our results mainly on Ha (on this sample we were able to 
reach larger waiting times). Hb has been used to check finite-size effects. 

Using SPt, we computed the 1RSB free energy density F(m, (3) and complexity £(T) = 
(3 d m F(m, (3)\ m= i for our samples as a function of the temperature T = 1/(3. The resulting 
complexity is reported in Fig. [HI for sample Ha- The dynamic and static temperatures are 
defined, respectively as the points where a non-trivial (1RSB) solution of the cavity equations 
first appears, and where its complexity vanishes. From the results of Fig. |S] (left frame) we 
get the estimates T d = 0.557(2) and T c = 0.467(2). 

In analogy with the analytic solution of the p-spin spherical model {241125] , we assume the 
aging dynamics of the model (|4.1|) to be dominated by threshold states. These are defined as 
the 1RSB metastable states with the highest free energy density. Although not exact we 
expect this assumption to be a good approximation for not-too-high values of 7. The threshold 
1RSB parameter mth(T) can be computed by imposing the condition c^[m.F(m, (3)] = 0. 
We computed m t h(T) on sample Ha for a few temperatures below T^. We get m t h(0.3) = 
0.395(10), m t h(0.4) = 0.58(1), m t h(0.5) = 0.80(1). Moreover, in the zero temperature limit, 
we obtain m t h(T) = HthT + 0(T 2 ), with fj, t ^ = 1.08(1). These results are summarized in 
Fig- El (right frame). A good description of the temperature dependence is obtained using the 
polynomial fit m t h(T) = 1.08T + 0.038T 2 + 2.17T 3 (cf. continuous line in Fig. |HJ right frame). 

Now we are in the position of precising the connection between single-spin statics and 
aging dynamics, outlined in Sec. El It is convenient to work with the integrated response 
functions Xi(t,tw)- Equation (|2.3j) implies the relation x«(Mu>) = Xi[Ci(t,t w )] to hold in the 
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limit t,t w — > oo. Within a 1RSB approximation, Eq. (|2.4I) corresponds to: 

^ r i J 1 - 9 for g > ^]. th , , s 

I, 1 — ?EA,th ~~ m th\Q — <?EA,thJ IOT Q < <?EA,th ! 

where we used the shorthand q^ th = (?EA( m th)- Since the SPt algorithm allows us to compute 
both rrith and the parameters q^lim) for a given sample in linear time, we can check the above 
prediction in our simulations. 



4.1 Numerical results 

We ran our simulations at three different temperatures (T = 0.3, 0.4, 0.5) and intensities of 
the external field (/io = 0.05,0.1,0.15). In order to probe the aging regime, we repeated 
our simulations for several waiting times t w = 10, 10 2 , 10 3 , 10 4 , with (respectively) AtMAX = 
2l3 2*6 2 16 2 18 

We summarize in the Table below the statistic of our simulations on sample TLa- 
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For sample He, we limited ourselves to the case ho = 0.10, T = 0.4, and generated 0.9 • 10' 
Metropolis trajectories with t w = 10 3 . 



4.1.1 Two types of spins 

The most evident feature of our numerical data, is that the spins can be clearly classified in 
two groups: (I) the ones which behave as if the system were in equilibrium: the corresponding 
correlation and response functions satisfy time-traslation invariance and FDT; (II) the out- 
of-equilibrium spins, whose correlation and response functions are non-homogeneous on long 
time scales and violate FDT. 

Of course the group (I) includes the isolated sites, but also an extensive fraction of non- 
isolated sites (for instance the 12 sites z = 1,6, 8, 14, 27, 39, 68, 74, 77, 84, 87, 98 of sample Ha)- 

(i) 

Remarkably these sites are the ones for which the SPt algorithm returns q^ A = 0: they are 
paramagnetic from the static point of view. 

In Figs. El and ^1 we present the correlation function and the FD plot, respectively, for a 
type-I site and a type-II site. In both cases we took T = 0.5 and h = 0.05. 

There exists a nice geometrical characterization of type-I sites in terms of a leaf-removal 
algorithm 41,42 . Let us recall here the definition of this procedure. The algorithm starts by 
removing all the interactions which involve at least one site with connectivity 1. The same 
operation is repeated recursively until no connectivity- 1 site is left. The reduced graph will 
contain either isolated sites, or sites which have connectivity greater than one. The sites of 
this last type are surely type-II, but they are not the only ones. In fact one has to restore a 
subset of the original interactions according to the following recursive rule. If an interaction 
involves at least two type-II sites, restore it and declare the third site to be type-II. If no such 
an interaction can be found among the original ones, stop. In this way, one has singled out the 
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1 10 100 1000 10000 100000 1e+06 0.2 0.4 0.6 0.8 1 

t-t w C 

Figure 9: Correlation function (left) and FD plot (right) of the spin % — 1 (sample Ha) for T = 0.5, 
ho — 0.1 and t w = 10 -rlO 4 . Time-traslation invariance is well verified for t w > 100. The discrepancy 
from FDT (continuous line on the right) can be ascribed to nonlinear response effects. 




Figure 10: Correlation function (left) and FD plot (right) for the spin i — (sample Ha)- The 
dashed line on the left corresponds the ergodicity breaking parameter <j4a = 0.716(7) obtained with 
the SPt algorithm. On the right we report with a full line the corresponding OFDR predicted 
within a 1RSB scenario. 
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Figure 11: The correlation functions of four type-I spins of sample H.^: from above to below i = 
1, 14, 8, 15. The first three sites are connected to the rest of the cluster (and therefore interacting), 
the last one is isolated (free). 



subset of original interactions which are relevant for aging dynamics. The sites which remain 
isolated after this restoration procedure are type-I sites, the others are type-II. 

The dynamical relevance of this construction is easily understood by considering two simple 
cases. A connectivity- 1 spin whose two neighbors evolve slowly will be affected by a slow local 
field and will relax on the same time-scale of the field. In the opposite case, two connectivity- 
1 spins whose neighbor evolve slowly will effectively see just a slowly alternating two-spin 
coupling between them. They will relax as fast as a two-spin isolated cluster does. 

It is worth stressing that the above construction does not contain all the dynamical infor- 
mation on the model. For instance, one may wonder whether the dynamics of type-I spins does 
resemble to the dynamics of isolated spins. The answer is given in Fig. ^Jwhere we reproduce 
the correlation functions for several different type-I spins, for T = 0.5, h = 0.1 and t w = 10 4 
(remember that the dependence upon t w is weak for these sites). The results are strongly 
site-dependent and by no way similar to the free-spin case. Notice the peculiar behavior of 
the isolated spin, an artifact of Metropolis algorithm with sequential updatings. Weren't for 
the perturbing field we would have Oi(t) = (— l)*<jj(0), which implies Ci(t w + l,t w ) = — 1, 
Ci(t w + 2, t w ) = 1, and Ci(t w + 2At, t w ) = for At > 2 [remind the time average in Eq. (|2.5|> ]. 
In the presence of an external field the correlation function (with no time-average) becomes 
Ci(t w + t, t w ) = [- exp(-2/?/i)]*. 

Finally, in Fig. ^]we compare the numerical results with the prediction from the statics, 
cf. Eq. (|4.2|) . The agreement is quite good although finite-i^ and finite-/io effects are not 
negligible. 

4.1.2 Dependence on the perturbing field 

We claimed that type-I spins satisfy FDT. However the FD plot in Fig. E] is not completely 
convincing on this point. The data show a small but significant discrepancy with respect to 
the FDT-line. We want to show here that this is (at least partially) a finite-/io effect. 

In order to have a careful check of this effect, we used, on sample TCa, three different values 
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Figure 12: The fluctuation-dissipation plot for a type-I (left frame: site i = 1 on Ha) and for a 
type-II site (right frame: site % = on Ha) at T = 0.4. Different symbols refer to different values of 
the perturbing field: h = 0.05 (x), 0.10 (*), 0.15 (□), and h — > extrapolation (continuous line). 
The dot-dashed lines refer to the FDT and (in the right frame) the predicted OFDR. 



of ho for any choice of T and t w . The main outcome of this analysis is that the ho dependence 
is much stronger for type-I than for type-II spins. This can be easily understood: type-II spins 
interact strongly with their neighbors and are therefore "stiffer". 

In Fig. ^]we present the FD plots for two sites of sample Ha'- i = 1 (type-I) and i = 
(type-II), and the three values of ho used in simulations. We expect the correlation and 
response functions to approach the zero-field limit as follows: 

Ci(t, t w \h ) = a t (t, t w ) + Di(t, t w ) hi + 0{ht) , (4.3) 
Xi(t, t w \h ) = Xi{t, t w ) + 5i(t, t w ) hi + OQii) . (4.4) 

We fitted our numerical data using the above expressions (and neglecting next-to-leading 
corrections). The results of this extrapolation are reported in Fig. Altogether with the 
finite-Zio data. 

As can be verified on Fig. I12( for type-II spins the ho — > limit is well approximated 
(within 2 -r- 3%) by the ho = 0.05, or (within 7 -r- 8%) even by the ho = 0.1 result. In what 
follows, we shall generally neglect finite-/io artifacts for these sites. On the other hand, an 
extrapolation of the type (|4.3|l - (|4.4|l is necessary for type-I spins. 

4.1.3 Glassy degrees of freedom 

In this Section we focus on type-II sites, which remain out-of-equilibrium on long time scales. 
In Fig. A] we reproduce the correlation and response functions of all the spins of sample H.a 
in a movie plot. We fix t w = 10 4 and watch the single-spin correlation and response functions, 
as the system evolves, i.e. as t grows. The behavior can be described as follows: (i) for small t, 
all the points (Ci,Xi) stay on the fluctuation dissipation line T\i = 1 — Q, type-I and type-II 
spins cannot be distinguished; (ii) as t grows, type-I spins reveal to be "faster" than type-II 
ones and move rapidly toward the C = 0,x = 1 corner; (Hi) just after this, type-II spins 
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' 0.0 0.5 1.0 0.5 1.0 

C(t,t ) 

Figure 13: Movie plot for sample Ha of Sec. HI Here we use h = 0.1, T = 0.4, and t w = 10 4 . The 
different frames correspond to (from left to right and top to bottom): At = 2 4 , 2 9 , 2 12 , 2 15 , 2 16 , 
and 2 17 . Black and white circles refer, respectively, to type-I and type-II sites. Continuous lines 
correspond to ordinary FDT x» = (1 — Cj)/T, while dotted ones are fits to a modified relation 
Xi = (1 - Ci)/T movie . We get T movie = 0.459 (for t = 2 12 ), 0.536 (t = 2 15 ), 0.564 (t = 2 16 ), 0.590 
(i = 2 17 ). 
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Figure 14: On the left, the FD plot for seven different sites: i = 10 (+), 11 (o), 12 (filled squares), 
13 (□), 16 (x), 17 (*), 18 (•). On the right, a collapse plot of the same data, cf. Eq. (|4.6j) . Here 
T = 0.3, h = 0.1, t w = 10 4 , and we use q = 0.94. 



move out of the FDT relation, all together 3 ; (iv) type-II keep evolving in the C-\ plane but, 
amazingly, they stay, at each time on a unique (moving) line passing through C = 1, % = 0. 
On the same graphs, in Fig. ^3 we show the results of a fit of the type 

Xi(t,t w ) = - 1 7 -—[1-C i (t,t w )}. (4.5) 

movie \"j ) 

The fit works quite well: it allows to define a new effective temperature: the "movie" temper- 
ature T mov i e (t,t w ). The thermometrical interpretation of T mov i e (t,t w ) will be discussed in Sec. 
El T movie (t,t w ) increases with t at fixed waiting time t w . Notice the difference between this 
formula, and Eq. (|.3.14[) which we argued to hold for coarsening systems. The organization of 
heterogeneous degrees of freedom in the \-C plane is strongly dependent upon the nature of 
the physical system as a whole. 

The cautious reader will notice a few discrepancies between the above description and the 
data in Fig. 1131 Type-I spins reach the (C = 0, x = 1) corner slightly after type-II ones move 
out of the FDT line. A careful check shows that this is a finite-^ artifact. Moreover, for large 
times, they stay slightly below the FDT line. This phenomenon has been discussed in Sec. 
I4.1.2l and proved to be a fmite-/io effect. 

Let us now consider the local OFDR's, and compare the dynamical results with the static 
1RSB prediction (|4.2|) . A preliminary check was given in Fig. 1101 In Fig. 1141 left frame, we 
reproduce the Xi versus Cj curves for 7 type-II sites. They are superimposed for short times 
(quasi-equilibrium regime) and spread at later times (aging regime), but remaining roughly 
parallel to each other. If the static prediction Q4.2J) holds, we can collapse the various x»[Ci] 
curves by properly rescaling x% an d Cj- A particular form of rescaling, which is quite natural 

3 Here expressions like "all together" must be understood as "on the same time scale in the aging limit". Let 
us, for instance, assume that two spins i and j begin to violate FDT after times U(t w ) and tj(t w ). We say that 
they begin to violate FDT "together" if lim^^oo h(ti)/h(tj) — 1 (h(t) being the appropriate time parametrization 
function). 
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Figure 15: Correlation between the theoretical prediction for the local EA parameters and the 
results of out-of-equilibrium simulations. On the left we show the data for sample Ha (T = 0.5, 
t w = 10 3 , and h = 0.1), on the right for sample H B (T = 0.4, t w = 10 3 , and h = 0.1). 
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Figure 16: Distribution of the slopes of single-site OFDR's for T = 0.4. The vertical lines correspond 
to the theoretical prediction for the 1RSB parameter m t h- On the left (right) we kept h = 0.05 
(t w = 10 4 ) fixed. 
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Figure 17: Movie plot for sample 1~Lb (T = 0.4, ho = 0.1): we show the position of all the degrees 
of freedom in the %-C plane, for t w = 10 3 and At = 2 16 . The thin continuous lines are the FD plots 
for a few selected sites (in this case At varies between and 2 16 ). In the inset: the histogram of 
slopes of the FD curves in the out-of-equilibrium regime. 



for coarsening systems, was used in Sec. 13.1.21 cf. Eqs. ()3.13j) . It turns out that, in this case, 
a better collapse can be obtained by using the definition below: 

cr = i- *~ ( f (i-Q), xr= m (4.6) 

1 ?EA,th 1 ?EA,th 

where q is a reference overlap (which can be chosen freely). In Fig. 1141 right frame, we plot 
xl es and Cf es for the same 7 spins as before, computing the c/ EAith with the SPt algorithm. 
Note that there is no fitting parameter in this scaling plot. 

It can be interesting to have a more general look at the statics-dynamics relation. In order 
to make a comparison, we fitted 4 the single-site Xi-versus-Cj data to the theoretical prediction 
(|4.2|1 . The results for the two fitting parameters and m^ i,&t \ are compared in Figs. 

I15| 1161 with the outcome of the SPt algorithm. Although several sources of error affect the 
determination of the EA parameters from dynamical data, the agreement is quite satisfying. 

In the above paragraphs we stressed two properties of the aging dynamics of the model 
(|4~Tj) : the alignment in the movie plots, cf. Fig. [T3]and Eq. ((OJ) . and the OFDR (|4~2~)) . Let 
us notice that these two properties are not compatible at all times (t,t w ). In fact we expect 
our model to verify the weak ergodicity-breaking condition lim^oo Ci(t, t w ) = 0. Therefore, 
in this limit, the alignment (|4.5|) cannot be verified unless the x% become site-independent. On 
the other hand, this would invalidate the OFDR (|4.2j) . 

One plausible way-out to this contradiction is that Eq. (|4,5|) breaks down at large enough 
times. How this may happen is well illustrated by the numerical data concerning sample TLb 
shown in Fig. 1171 It is quite clear that the simple law (|4.5|) no longer holds. Nevertheless, 

4 While fitting the data we restricted ourselves to the C'i(t + t Wl t w ) < q^l time range. In practice we took t > 2 7 , 
2 10 , and 2 13 , respectively for t w = 10 2 , 10 3 , and 10 4 . 
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Figure 18: Qualitative picture of aging dynamics for discontinuous glasses. The black circles corre- 
spond to three different spins for a given pair of (large) times (t,t w ). Notice that, for intermediate 
times At ~ t w they stay on the same line passing through the point (C = l,x = 0). As t grows, 
they move with parallel velocities (arrows). Along the time, each of them describes a different 
fluctuation-dissipation curve (dashed lines). 



it remains a very good approximation for the sites with a large EA parameter q^l > 0.5. 
Moreover, it seems that the points corresponding to different sites still lie on the same curve 
in the x~C plane, although this curve is not a straight line as in Eq. (|4.5|) . We shall further 
comment on this point in Sec. [7J 

The general picture which holds at intermediate times (or large (?ea's) for discontinuous 
glasses is summarized in Fig. This should be compared with Fig. ^ which refers to 

coarsening systems. 



The Viana-Bray model ^T7\ is a prototypical example of continuous spin glass. It is defined by 
the Hamiltonian 



where the graph Q is constructed by randomly choosing M among the N(N — l)/2 couples 
of spins, and the couplings Jij are independent identically distributed random variables. The 
average connectivity of the graph is given by c = 2M/(N — 1). If we assume that the couplings 
distribution is even, the phase diagram of this model is quite simple [171143] . For c < 1 the 
interaction graph does not percolate and the model stays in its paramagnetic phase at all 
finite temperatures. For c > 1 the graph percolates and the giant component undergoes a 
paramagnetic-spin glass phase transition. The critical temperature is given by the solution of 
the equation Ej(tanh/?J) 2 = 1/c. Below the critical temperature, a finite Edwards- Anderson 
parameter q EA develops continuously from zero. 



5 Continuous glasses 




(5.1) 
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We considered three samples of this model: hereafter they will be denoted as Ga, Gb and 
Gc- The interaction graph and the signs of the interactions Jj,- were the same for Ga and 
Gb'- m particular we used N = 1000 and M = 1999, i.e. c ~ 4, and chosen the interaction 
signs to be ±1 with equal probabilities. The two samples differ only in the strength of the 
couplings. While in Ga we used \ Jij\ = 1, in Gb we took \ Jij\ = kJo, where k £ {1, . . . , 10} 
with uniform probability distribution and Jo = 0.161164 5 . We made this choice in order to 
check the effects of degenerate coupling strengths on the aging dynamics. The sample Gc was 
instead much larger: we used ./V = 10000, M = 20190 (once again c ~ 4) and Jij = ±1 with 
equal probabilities. The critical temperatures for c = 4 and the two coupling distributions 
used here are T c m 1.8204789 (for Ga and Gc) and 1.6717415 (Gb)- 

The glassy phase of the VB model is thought to be characterized by FRSB. Nevertheless we 
can use the SPt algorithm to compute a one-step approximation to the local overlaps and the 
local OFDR's. Of course, such an approximation will have the simple two time-sectors form, 
see Eq. (|4.2|) . instead of the expected infinite-time-sectors behavior. However the situation is 
not that simple because of two problems: 

• We expect, in analogy with the Sherrington-Kirkpatrick model f2H|, the dynamics of this 
model to reach the equilibrium free-energy in the long time limit. It is not clear whether 
a better approximation to the correct OFDR is obtained by using the threshold value 
m t h or the ground-state value m gs of the 1RSB parameter. 

• The SPt algorithm does not converge. After a fast transient the probability distributions 
of local fields oscillate indefinitely. This is, plausibly, a trace of FRSB. 

The first problem does not cause great trouble because the two determinations of m are, 
generally speaking, quite close. On the other hand, we elaborated two different way-out to the 
second one: (i) To force the local field distributions to be symmetric (which can be expected to 
be true on physical grounds): this assures convergence, (ii) To average the local EA parameters 
over sufficiently many iterations of the algorithm. 

While the approach (i) seems physically more sound, it underestimates grossly the Q^^'s. 
The approach (ii), which will be adopted in our analysis, gives much more reasonable results. 
Notice that the authors of Refs. |22[ I44| followed the same route. In their calculation, they 
faced no problem of convergence. In fact they required convergence in distribution, while we 
require convergence site-by-site. 

5.1 Numerical results 

Most of our simulations were run at temperature T = 0.5, and with /in = 0.1. We used waiting 
times t w = 10 2 , 10 3 , 10 4 , and, respectively At M AX = 2 14 , 2 16 , 2 18 . In the Table below we 
summarize the statistics used in each case. 
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5 This distribution has a three nice properties: (i) being discrete, it allows to precompute the table of Metropolis 
acceptance probabilities; (ii) since we will use temperatures T 3> Jo, it is a good approximation to a continuum 
distribution; (Hi) it has Ej(J 2 ) w 1. 
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Figure 19: Single-spin correlation and response functions for the sample Qa (VB model) for T = 0.5, 
ho = 0.1, and t w = 10 4 . The continuous line and full circle refer to the global correlation and 
response. 



Moreover we simulated iV s t a t = 4.2 • 10 5 Metropolis trajectories at temperature T = 0.4 on 
sample Qa with t w = 10 4 and A^max = 2 18 . 

In Fig. EI]we show the movie plot of sample Qa for t w = 10 4 . As in the previous Sections, 
the local correlation and response functions are strongly heterogeneous: the global two-times 
functions give just a rough idea of the dynamics of the system. Moreover all the points quit the 
FDT line on the same time scale in the aging limit (cf. Sec. 14.1.3(1 . However, their behavior 
in the aging regime does not fit any of the alignment patterns we singled-out in the case of 
coarsening systems, cf. Eq. (|3,14|) and Fig. ^ or discontinuous glasses, cf. Eq. (|4.5|) and Fig. 
ITSl We repeated the same type of analysis for the numerical data obtained on sample Qc- 
In this case, see Fig. 1201 the points corresponding to local correlation and response functions 
are much less spread in the \~C plane. Therefore our simulations are quite inconclusive on 
the possibility of defining a "movie" temperature as in Eq. (|4.5j) . To settle the question, 
simulations on larger samples are probably necessary. 

Numerical results on sample Qa are also deceiving for what concerns local OFDR's, cf. Fig. 
1211 It seems that the local FD plots depend strongly upon the waiting time and the particular 
site. Moreover the slopes of this plots (for a given couple (t w ,At)) change from site to site. 
These effects are much smaller in sample Qc ■ In Fig. 22 we consider the distribution of slopes 
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Figure 20: Single-spin correlation and response functions: here we compare the results obtained on 
samples Qa (°) and Q c (•) which are of different sizes. 




Ci(t,t w ) Ci(t,t w ) 

Figure 21: FD plot for a few selected sites of sample Qa (T = 0.5, ho = 0.1). Notice the completely 
different behaviors of the sites in the two frames. The sites in the left frame, with connectivity 1 
(site 111) and 3 (site 164), look like a "glassy" system. The ones on the right, with connectivity 4 
(site 103) and 2 (site 114), look like a "coarsening" system. 
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Figure 22: Distributions of slopes of the local FD plots in the aging regime. The two curves refer 
to sample Ga (°) an d Gc (•)> which are of different sizes. 



of local FD plots for samples Ga and Gc- We computed the slopes by fitting the aging part of 
the plot to the one-step form (|4.2|) . 

By the same fitting procedure we extracted the local EA parameters. The comparison 
with the predictions of the SPt algorithm, cf. Fig. |23]is quite satisfying. Notice that, both 
in analyzing the numerical data and in using the SP^ algorithm, we are adopting a 1RSB 
approximation, cf. Eq. (|4.2|) . to the real OFDR. The slopes considered in Fig. I22I should 
therefore be understood as average slopes in the aging regime. We expect the systematic error 
induced by this approximation to be small. 

The arguments of Ref. [13] imply that the slopes (effective temperatures) of the OFDR's 
for different degrees of freedom should be identical. This conclusion is valid only in the aging 
window 1 <C t w ,At <C t erg (N). Our numerical data, cf. Fig. E21 suggest a clear trend 
confirming this expectation. Nevertheless, they show large finite-size effects due, arguably, 
to a mild divergence of t erg (iV) with N: the smaller (N = 10 3 ) samples begin to equilibrate 
during the simulations. This is quite different from what happens with discontinuous glasses, 
cf. Sec. H In that case, we did not detect any evidence of equilibration even in sample TCa 
(N = 10 2 ). A better understanding of the scaling of t erg (N) in different classes of models 
would be welcome. 

6 Weakly interacting spins 

We lack of analytical tools for studying the dynamics of diluted mean-field spin glasses (for 
some recent work see Refs. |46U47lH51l49j ) . This makes somehow ambiguous the interpretation 
of many numerical results. For instance, the identity of effective temperatures for different 
spins, although consistent with our data, see Figs. EDandEU could still be questioned. This 
would contradict the general arguments of Refs. |45|I5()|. Even more puzzling is the definition 
of a "movie" temperature along the lines of Eq. (|4.5|) . Such a definition seems to be consistent 
only in some particular models and time-regimes. In this Section we want to point out a simple 
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Figure 23: Local EA parameters on sample Qc- Numerical results, obtained by fitting the aging 
part of the FD plots, are compared with the outcome of the SP^ algorithm. 



perturbative calculation which supports the identity of single-spin effective temperatures, in 
agreement with the standard wisdom. Moreover it give some intuition on the range of validity 
of the definition (|4.5j) . 

Let us consider a generic diluted mean-field spin glass with fe-spin interactions: 

H (°) = - J « (Jq i ' • • • ' °ak • (i 6 - 1 ) 

aen 

Here a = {ai, . . . is a fc-uple of interacting spins, and H is a fc-hypergraph, i.e. a set of 
M such fc-uples. 

Let us focus on a particular site, for instance i = 0, and assume that it is weakly coupled 
to its neighbors. It is quite natural to think that its response and correlation functions can 
be related to the response and correlation functions of the neighbors. To the lowest order this 
relation reads: 

Co g (t,t w ) = ]T(tanh/3J„) 2 \\ C* s (t,t w ) + 0([3 4 J 4 ) , (6.2) 

Q30 iGa\0 

RoHt,t w ) = 5>anh/3J„) 2 Rt g (t,t w ) J] C; e (t,t w )+0((3 4 J 4 ). (6.3) 

a30 i£a\0 jG«\{CM} 

We shall not give here the details of the derivation. The basic idea is to use an appropriate 
dynamic generalization of the cavity method |511l52j . As for static calculations |21| . this 
approach gives access to single-site quantities for a given disorder realization. Notice that Eq. 
(|6.2[) can be easily obtained by assuming that the spin <jq does not react on its neighbors. This 
is not the case for Eq. (|6.3|) , 

Equation (j6.2|) implies a relation between local Edwards- Anderson parameters: 

?S = ^(tanh/3J„) 2 J] q V +0 (p 4 J 4 ). (6.4) 

Q30 i€a\0 
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Figure 24: FD plots for a few weakly-interacting spins: numerical results (•) and outcomes of 
the perturbative formulae (J6.2j) and (J6.5j) (o). We consider connectivity- 1 sites on the left, and 
connectivity-2 sites on the right. 



In the k = 2 (Viana-Bray) case, we can derive from Eq. (|6.3|) a simple relation between the 
integrated responses: 

l-T X f(t,t w ) = ^(tanh/?J 0i ) 2 [l-r X f(t,^)]+O(/3 4 J 4 ), (6.5) 

iedO 

where dio denote the set of neighbors of the spin io- In the general (k > 2) case Eq. (j6.3[) 
cannot be integrated without further assumptions. 

We checked the above relations on our numerical data for the Viana-Bray model. Sample 
Qb is particularly suited for this task, since we can choose spins whose interactions have a 
varying strength. In Fig. 1241 we consider a few spins with connectivity one and two, and 
compare their correlation and response functions with the outcome of Eqs. (|6.2j) and 1)6. 5 J) . Of 
course, the perturbative formulae are well verified only for small couplings. For connectivity-2 
sites we have plotted in Fig. only those with coupling of the same strenght, since spins with 
2 couplings of very different strengths behave very similarly to connectivity- 1 spins. 

Let us now discuss some implications of Eqs. 1)6.2)1 . (|6.3|) . If we define the fluctuation- 
dissipation ratio as Xi(t,t w ) = TR^ g (t,t w )/dt w C^ g (t,t w ), we get: 

X (t, t w ) « = , (6.6) 

where 

Wa,i(t,t w ) = (tanh /3Ja) 2 d tw Ci(t,t w ) J] Cj(t,t w ) (6.7) 

iea\{i,o} 

are positive weights. Therefore, at the lowest order in perturbation theory, the effective tem- 
perature of the spin o"o is a weighted average of the effective temperatures of its neighbors. 
Let us suppose that this conclusion remains qualitatively true beyond perturbation theory. It 
follows that Xi(t,t w ) = X(t,t w ) is independent of the site i. In fact, if the Xi(t,t w ) were 
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site-dependent we could just consider a site such that Xi t (t,t w ) is a relative maximum and 
show that Eq. ([6.6)) cannot hold on such a site. With a suggestive rephrasing we may say that 
effective temperatures must diffuse until they becomes site- independent. 

Moreover, Eqs. ([6.2)1 and (|6.3j) can be used to construct examples of weakly interacting 
spins which violate the alignment in the x~C plane which we encountered for discontinuous 
glasses, cf. Eq. ([4.5)1 and Fig. 1131 The simplest of such examples is obtained by considering 
the Viana-Bray (A: = 2) case, and assuming that the site has just one neighbor. In this case 
it is immediate to show that 

C g (t, t w ) ^ C i g (t, t w ) 
l-T Xo (t,t w ) ~ l-T Xi (t,t w ) ' { ' ' 

i.e. weakly interacting spins have the tendency to align as in coarsening systems. The reader 
can easily construct analogous examples for k > 2 models. This suggests that the "movie" 
temperature (|4.5|) is well defined uniquely for strongly interacting and glassy systems, or, in 

(i) 

other words, for slow-evolving sites with a q^i close to one. 



7 Discussion 

In the last two Sections we shall discuss the general properties of single-spin correlation and 
response functions which emerge from the numerics. We will give an overview of such properties 
in the present Section, and precise the thermometric interpretation of some of them in the next 
one. 

We shall focus on two-times correlation and response functions Ci(t,t w ) and Ri(t,t w ) (see 
[4 7) for a preliminary discussion of multi-time functions) and distinguish two types of facts: (i) 
their scaling behavior in the large time limit; (ii) the fluctuation-dissipation relations which 
connect correlation and response. 

7.1 Time scaling 

Following Refs. [2^1130], we assume monotonicity of the two-times functions: 3jCj(i, t w ), 
dtRi(t,t w ) < 0, and dt w Ci{t,t w ), dt w Ri(t,t w ) > 0. Moreover we consider a weak-ergodicity 
breaking situation: Ci(t,t w ), Ri(t,t w ) — > as t — ► oo for any fixed t w . All these properties are 
well realized within our models. 

It is quite natural to assume 6 that, for pair of sites i and j there exist two continuous 
functions and fji such that 

Ci(t,t w ~) — fyj [Cj (t , iui)] > Cj(t 1 t VJ ) — fji [Ci (t, ito)] , 

in the t,t w — > oo limit. Notice that we can always write 

Ci(t,t w ) = fij[Cj(t,t w ),t] . (7.2) 

We are therefore assuming that the functions [C, t] admit a limit as t — > oo and that the 
limit is continuous. Since fij[C, t] is smooth and dcfij[C,t] > 0, if the limit exists it must be 

6 A similar statement appears in Ref. |55j . Notice however two differences: the cited authors consider correlation 
functions for different modes in Fourier space, while we consider different sites; more important, they take the average 
over quenched disorder, while we consider a fixed realization of the disorder (our statement would be trivial if we 
took the disorder average). 



(7.1) 



29 



a continuous, non-decreasing function of C. Since Eq. (|7.1[) implies that both and fji are 
invertible (indeed fij o = 1, see below) they must be strictly increasing. 

Without any further specification, the property l|7.1|) is trivially false. Consider the example 
of type-I (paramagnetic) spins in the 3-spin model studied in Sec. If i is type-I, and j is 
type-II, Ci(t,t w ) — > in the aging regime, while Cj(t,t w ) remain non-trivial: /«[•] cannot be 
inverted. Another example, would be that of a Viana-Bray model, cf. Sec. 03 such that the 
interactions graph has two disconnected components. 

However, both these counter-examples are somehow "pathological". We can precise this 
intuition by noticing that Eq. (|7.1|) defines an equivalence relation (in mathematical sense) 
between the sites i and j. Therefore the physical system breaks up into dynamically connected 
components which are the equivalence classes of this relation. Type-I and type-II spins in the 
3-spin model of Sec. |l]are two examples of dynamically connected components. Hereafter we 
shall restrict our attention to a single dynamically connected component. Physically, structural 
rearrangements occurs coherently within such a component. 

Clearly the transition functions {fij} have the following two properties: (£) fji = f^j , and 
(**) fij = fik ° fkj- This implies that they can be written in the form fij = ff o fj (the proof 
consists in taking a reference spin k = and writing fij = fa o f j = f~ x o / 0j ). Of course the 
functions fi are not unique: in particular they can be modified by a global reparametrization 
fi-> 9° fi- 

Although very simple, the hypothesis (|7.1j) has some important consequences. Suppose 
that Cj(t,t w ) has p discrete correlation scales (in the sense of Refs. |29| l30]^. characterized by 
Qa+l < Cj(t, t w ) < qi\ for a = 1, . . . ,p. Within a scale we have 

C 3 {t,t w )^cf\h^(t)/h^(t w )). (7.3) 

where ha (t) is a monotonously increasing time-scaling function. Two times t and t w belong 
to the same time sector if 1 < ha (t)/ni (t w ) < oo. 

Applying the transition function fij to the above equation, one can prove that, for each 

scale a of the site j, there exists a correlation scale for the site i, with q^ +l < Ci(t,t w ) < q$ 

and q. = /«[?. ]. Moreover ha (t) = ha (t) = h a (t) (up to an irrelevant multiplicative 
constant) and 

Cl a) = fij o Cf . (7.4) 

In summary there is a one-to-one correspondence between the correlation scales of any two 
sites. Notice that this is a necessary hypothesis if we want the connection between statics and 
dynamics 27 , 28 to be satisfied both at the level of global and local (single-spin) observables. 
A spectacular demonstration of the correspondence of correlation scales on different sites is 
given by our movie plots, cf. Figs. 03 El and 1191 In particular such correspondence implies 
that all the (xi, Q) points leave the FDT line at once. 

Equation ()7.1|) can be rephrased by saying that the behavior of one spin "determines" the 
behavior of the whole system. This is compatible with the locality of the underlying dynamics 
because: {%) "determines" has to be understood in average sense; (u) the relation ()7.1|) is not 
true but in the aging limit. 

7.2 Fluctuation-dissipation relations 

On general grounds, we expect single-spin quantities to to satisfy site-dependent OFDR's of 
the type (|2.3|) . In integrated form we obtain, for large times t,t w , the relation Xi(Mw) = 
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Xi[Ci(t, t w )\. We think that we accumulated convincing numerical evidence in this direction as 
far as the models of Sees. (coarsening) and E] (discontinuous spin glass) are considered. The 
situation is more ambiguous (and probably very hard to settle numerically) for the Viana-Bray 
model of Sec. 03 

Fluctuation-dissipation relations on different sites are not unrelated: we expect [IS] to be 
able to define a site-independent effective temperature as follows 



numerics support this identity both for coarsening systems, cf. Sec. 03 and discontinuous 
glasses, cf. Sec. For continuous glasses, cf. Sec. 03 the situation is less definite. In Sec. 03 
we presented a perturbative calculation which support Eq. (|7,5|) also in this case. 

A suggestive approach [SU] for justifying (|7.5|) consists in regarding T e fj(i, t w ) as the tem- 
perature measured by a thermometer coupled to a particular observable of the system. It is 
quite natural to think that the result of this measure should not depend upon the observable. 
In aging systems with more than just one time sector, this approach is not consistent unless 
the following identity holds: 



Proving this conclusion will be the object of the next Section. The new effective temperature 
7movic(^ t w ) is in fact the one measured by a particular class of thermometers which we shall 
denote as "sharp". It is a weighted average of the effective temperatures (in the sense of <|7.5JI ) 
corresponding to different time sectors. 

Let us notice that Eq. (|7.6[) is remarkably well verified in our discontinuous spin-glass 
model, cf. Fig. although it breaks down for (t/t w ) ^> 1. In Sec. 01 we demonstrated that 
it does not hold for coarsening systems, and in fact a different relation is true in this case, 
cf. Eq. (|3.14|) . Finally, we were not able to reach any definite conclusion for the Viana-Bray 
model of Sec. 03 

8 Thermometric interpretation 

One of the most striking (and unexpected) observations we made on our numerical results is 
summarized in Eq. (|7.6|) . In this Section we shall try to connect this empiric observation to 
thermometric considerations. This type of approach was pioneered in Ref. j^O] with the aim 
of interpreting the OFDR in terms of effective temperatures. Here we will show that, in a 
system with more than a single time scale, this interpretation is well founded if and only if the 
"movie" temperature ()7.6|) can be defined. In order to obtain this result, we shall carefully 
reconsider the arguments of Refs. [50 1 153 1 153] ■ 

According to Ref. jHH], the temperature of a out-of-equilibrium system can be measured 
by weakly coupling it to a "thermometer" , i.e. to a physical device which can be equilibrated 
at a tunable temperature Tth = 1 /Pth ■ The temperature of the system is defined as the value 
of T t h such that the heat flow between it and the thermometer vanishes. The details of the 
thermometer are immaterial in the weak-coupling limit. What matters are the correlation and 





(7.6) 
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response functions of the thermometer 7 C^(t, t w ) = C t h{t — t w ) and Rth{t,t w ) = Rth(t — t w ), 
which are assumed to satisfy FDT: i? t h( r ) = — Ah<9-7-Cth( T )- 

In the spirit of our work, we shall couple the thermometer to a single spin variable 
between times and t, and average over many thermal histories. The measured temperature 
Pth is given by 

th [ dt w R th {t - t w )d tw C t (t,t w ) = [ dt w R th {t - U(-Xi[^(M«,)])3 tl „Q(t,U , (8.1) 
J o J o 

where we assumed the general OFDR (|2,3|) in its integrated form: Xi{t-> tw) = Xi[Ci(t, t w )], and 
denoted by a prime the derivative of Xi['] with respect to its argument. Notice that a priori 
the measured temperature depends upon i and t, for a given thermometer. 

It is convenient to change variables from t w to q = Ci(t,t w ). This relation can be inverted 
by defining the time scale Tj(i; q) as follows: 

Ci{t,t-T t (t;q)) =q. (8.2) 

Using these definitions in Eq. (|8.1|) . we get 

Ah/ dqR th (r i (t;q))= [ dq R th (n(t; q^-xM) , (8-3) 



nun 



where g m i n = Ci(t, 0). As t — > oo, we have g m i n — > 0. In the same limit Ti(t;q) — > T 4 eq (g) if 

9 > ?ba, while g) -> oo if g < q® A . 

In order to measure temperatures on long time scales, we need a thermometer with an 
adjustable time scale. Mathematically speaking, we take Rth(T~) = Rth{^ I Vth) ? an d use r t h to 
select the time scale. The precise form of Rt^x) is not very important. We shall assume that 
Rth(x) ~ 1 for x <^ 1 and i? t h(x) ~ for x 3> 1. A simple example is Rth(%) = 0(x)e~ x . 
Some of the relations we will derive simplify if R t ^(x) ~ 9(x)9(x if — x). We will call such a 
thermometer "sharp". 

We have two type of choices for the thermometer time scale r^: 

1. We may take a "fast" thermometer, whose relaxation is much faster than the structural 
rearrangements of the system. Equivalently: we look at our thermometer after a time 
t ^> Tth- Mathematically this corresponds to taking the limit t — > oo with rth fixed. The 
result of such a measure is (for large times t) the bath temperature. 

2. We may use a "slow" thermometer, with a relaxation time which is of the same order of 
the time needed for a structural change in the system. This corresponds to taking the 
limits t — > oo, r t h ^ oo at the same time. If the system ages, the outcome of such a 
measure will depend upon the precise way these limits are taken. 

Let us consider separately the two cases. 

8.1 "Fast" thermometer 

In this case we have, as t — > oo, 

RMt; q)) - ( Fl{q) = ^(^WAth) for q > g|, (jU) 

for q < q^{, 



7 More precisely: Cth(t,t w ) and Rth{t,t w ) are the correlation and response functions for the observable of the 
thermometer which is coupled to the system. 
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with Fi(q^ A ) = and F;(l) = 1. Inserting into Eq. (jH3J we get 

Ah [ r dqF i (q)= [ dqF i {q)(-x' i [q))- (8-5) 

Assuming that in the "quasi-equilibrium" time sector (i.e. for Ci(t,t w ) > q^ A ) the system 
satisfies FDT, we can use x'i[o\ = ~Pi which yields /? t h = P-, as expected. 

8.2 "Slow" thermometer 

Here we shall assume that the system has p discrete correlation scales in the aging regime |3Uj . 
The generalization to a continuous set of correlation scales is straightforward. To each scale 
a € we associates a time-scaling function h a (t). As discussed in Sec. [7J h a (t) is 

site-independent. 

In order to probe the correlation scale a, we tune the thermometer time scale with the 
function Tth,A(t)- This function is defined by imposing 

ka ® -A, (8.6) 



h a (t - T thyA (t)) 

for some fixed number A > 1. 

Within the scale a, we have q^\_i < Ci(t, t') < qi 9 . It is easy to show that 

lim R th (n(t; q)/T thA (t)) = F i>a (q) , (8.7) 

with F i>a (q) = for q < q^ +1 , F i)0l (q) = 1 for q > q$ and F i>a (q) increasing in )■ 
Integrating by parts Eq. (|8.3|) . we get 

P$ I dqFl a (q)(l - q) = f dq F^q) Xi (q) , (8.8) 
Jo Jo 

which is our final expression for the temperature measured on the spin i (here we emphasized 
the dependence of /3th upon the site). 

Notice that the support of F[ a {q) is contained in the interval (ql^ +1 , q^). The expression 

(|8.8|) simplifies in two cases: (i) if the a-th correlation scale is "small" q^ +l ~ g4 ~ q* (and, 
in particular, when there is a continuous set of scales); (ii) if the thermometer is "sharp" in 
the sense defined at the beginning of this Section, and, therefore, F( a (q) is strongly peaked 

around some q*\ In both cases we have 

1 - q* 

Let us now imagine to couple two copies of the same thermometer to two different sites % and 
j. We shall measure two temperatures P$ « Xi(<7*^)/(1 ~ an d Pth ~ Xj(?* )/(l ~ I***), 
with g* = fij(qi )■ These two temperatures coincide, P^ sa p^ , only i/Eq. (J7.6|) is satisfied. 

The conclusion of the arguments presented so far is that the condition (|7.6j) is necessary 
if we want a given thermometer to measure the same temperature on any two spins of the 
system. Moreover this condition is sufficient for the special class of "sharp" thermometers. In 
the last part of this Section we will show that the condition Q7.6B is indeed sufficient for any 
thermometer, once ()7.5|) is assumed. 
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8.3 Thermometric equivalence of different sites 

We want to prove that Eqs. (|7.6|) and ()7.5j) imply the identity of thermometric temperatures 
on the sites i and j for any given thermometer. Let us stress that the measured temperature 
may, eventually, depend upon the thermometer. The essential ingredient for the "small intropy 
production" scenario of Ref. L 56 to be applicable, is that the result should not depend upon 
the site. 

Notice that from the definition ()8.2I) . it follows that the time scales defined on different 
sites are related as follows 

n(t;f ij (q)) = T j (t;q), (8.10) 

whence we easily derive the identity Fi t0l (fij(q)) = Fj^ a (q). By the change of variables q — > 
fij(q) we get, from Eq. (|8.8|) 

. q (J) -girt 
P® /J dqF' ha {q){1-nM))= /J dqF< ;a (q)Xi(fiM)> (8-H) 



a + l 



where we specified the range of q such that F'- a (q) is (possibly) nonzero. If use Eq. ()7.6|) to 
connect the responses on different sites, we obtain 



CO 



Ah dqF' ha (q)(l-q) 



U) 

a + l 



CO 



;;) dq F' j a {q)xj{q) 

a + l 



1 - fij(q) 



.12) 



The factors (1 — fij(q))/(l — q) prevent us from concluding that f3$ = [3$ with no further 

assumption. Let us assume Eq. (J73J), and that x'iW\ stays constant for q^ +1 < q < q^a ■ It 

follows that, within the scale a, fij(q) = 1 — /° ^ (1 — q), f®^ being a constant. This implies 

(i) (i) 
/? th = /3 t j^ for any thermometer. 
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A Large-n calculations 

In this Appendix we sketch the large-ra calculations whose results were presented in Sec. 13.1.11 

A.l Statics 

The trick for solving the periodic model of Sec. 13.1.11 is quite standard |58j . We define the 
n ■ V-components vector if) x which contains n components for each type of spin: 

r x ' u = €-i +u , ueA, xez\ (A.l) 

where x ■ I = X^=i an d A is the elementary cell. In this basis the Hamiltonian reads 

X,fl X 
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where K 



(p) 

au,bv 



bL u „ and 



The equilibrium correlation functions are computed by standard methods: 



where the V x V matrix Wl*(p) is given by 

d 



+ C* 



(A.3) 
(A.4) 



(A.5) 



The V Lagrange multipliers and the V magnetizations M u must be computed from the set 
of 2V equations given below: 



1 



£C(0)M e = 0, 

BZ( 2 ^) 



(A.6) 

(A.7) 



These equations have two type of solutions: at high temperature M u = and the matrix M^O) 
has rank V; at low temperature M u > and the matrix M*(0) has one vanishing eigenvalue. 

In the following Section we shall treat the dynamics of this model. Remarkably all the 
complication produced by inhomogeneous couplings affects the aging dynamics only through 
the values of the local magnetizations {M u }, the critical temperature T c and one more constant, 
A, which we are going to define. Consider the lowest lying eigenvalue A°(p) of the matrix M*(p). 
As p — > the corresponding eigenvector coincide with M v and A°(p) — ► 0. We then define 



A = Det 



~d 2 \°( P ) 




dp^dp v 


p=0_ 



(Ai 



All these quantities can be easily computed once the solution to Eqs. (|A.6D . (|A.7|) is known. 



A. 2 Dynamics 



The Langevin equation ()3.7|) is easily solved by defining the new order parameter rj) x as in the 
previous Section, going to Fourier space: 



(A.9) 



veA 



The "mass" matrix M™(p, t) is given by the expression (|A.5|) with the Lagrange multipliers 
replaced by their time-dependent version ( u (t). Of course lim^oo ( u (t) = Q . 
The correlation and response functions for the field tp x become V x V matrices. Their 

diagonal elements are the on-site correlation and response functions of the field 4>. Standard 

manipulations yield: 



C(t, t v 
R{t, t v 



d d p 

bzW? 
d d p 



U(p;t) 



tw 

1 + 2T / dsU{p]s)- l U{-p\s) 
o 



-i 



BZ 



(2tt) £ 



U(p;t)U(p;t) 



-i 



U{-p-t) (A.10) 
(A.ll) 
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The matrix U(p; t) satisfy the differential equation 

d t U(p; t) = -M(p; t)U(p; t) , U(p; 0) = I , 

and the Lagrange multipliers must satisfy the self-consistency conditions C uu (t,t) = 1. 
One can find the following asymptotic behavior for U(p;t): 



U(p; t) = At d '\l + -ft- d/2+1 + ...)e 
The constants A and 7 are simple numbers given below: 



U(p)t 



A 



7 



E M gA M u fq W4 A l/4 

1 + r/r* v ; 



r(i-d/2) 2 



(E„ eA M u 2 )(8vr)^AV2 r(2-d) 



The constant T* appearing in Eq. (jA. 14f) is defined as follows 

1 

— = J^dt &-u(0;t)- 2 a, 



(A.12) 

(A.13) 

(A.14) 

(A.15) 
(A.16) 

(A.17) 



where <r is the ^-dimensional unit vector parallel to the vector of the magnetizations: a u = 
M u / i^2 u M 2 ) 1 / 2 . The expression (|A.17|) is quite hard to evaluate, but this is not a problem, 
because T* cancels out in all physical quantities. 

Using the results listed above one can recover the general form (|3,3j) and the expressions 
(|3.1U|) - (|3.12|) . The universal functions which determine the domain wall contributions are given 
below for general dimension 2 < d < 4 (we recall that in the n — > 00 limit the model is well 
defined in non-integer dimensions): 



^c(A;d) 



r(i-d/2) 2 fx + x 
r(2 - d) 



>-i 



-d/2 



[1 + x 2 - d ) 



TJX;d) 



_ x d/2 



i-d+2 



1 + A 5 



2\-d+l / .2(1+A 2 )- 1 



dxx- d / 2 {l-x)- d l 2 , 



dxx d ' z (x d i 2 - \){x 2 - \y d i 2 . 



(A.18) 
(A.19) 



The integral in Eq. (jA.18|) diverges for d > 2: it is understood that it has to be analytically 
continued [SU] from d < 2 to obtain the correct result. 

It can be useful to consider the asymptotic behavior of the expressions (|A,18j) and l|A.19|) . 
As A — ► 00 (i.e.i 3> t w ) we have 



F C (X;d) 
T x (X;d) 



nd/2-1 \-d/2 
x -d+2 



T(l-d/2) 2 ^ 4 



r(2 - d) d-2 

r(i - d/2)r(i - d/4) 



+ r L" d/ ? 2 A-^ + o(A-'; 



+ 



1 



d-2 



r(2 - d) 
2 



(A.20) 



4-d 



2r(2 - 3d/4) 

As already remarked in Sec. l3.1.1l both functions vanish in the A — > 00 limit. 



A -2+d/2 + 0(A -3 + d/2 ) | 



(A.21) 
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When A — > 1 one gets 



^c(A;d) 



d 



(a - \y d / 2+1 



[1 + 0(A-1)] 



(A.22) 



2-d 



F x {\;d) 



2- d / 2 d 
4-d 



( A _!)2-d/2 [i + ( A _i)] . 



(A.23) 
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